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Abstract. We present a comprehensive asteroseismic study of the double-mode high-amplitude 6 Scuti star HD 187642 (RV 
Ari). The modelling includes some of the most recent techniques: 1) effects of rotation on both equilibrium models and adiabatic 
oscillation spectrum, 2) non-adiabatic study of radial and non-radial modes, 3) relationship between the fundamental radial 
mode and the first overtone in the framework of Petersen diagrams. The analysis reveals that two of the observed frequencies 
are very probably identified as the fundamental and first overtone radial modes. Analysis of the colour index variations, together 
with theoretical non-adiabatic calculations, points to models in the range of [7065,7245] K in effective temperature and of [1 190, 
1270] Myr in stellar age. These values were found to be compatible with those obtained using the three other asteroseismic 
techniques. 
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1. Introduction 

The 6 Scuti stars are intermediate-mass pulsating variables. 
With spectral types ranging from A2 to F0, they are located on 
and just off the main sequence in the lower part of the Cepheid 
instability strip (luminosity classes V & IV). Nowadays, the 6 
Scuti stars are considered as particularly suitable to the astero- 
seismic study of poorly known hydrodynamical process occur- 
ring in stellar interiors, such as convective overshoot, mixing 
of chemical elements, and redistribution of angular momentum 
1992), to mention only the most important. Due to the 
complexity of the oscillation spe ctra, their pulsating behaviour 
is not understood very well: see lCoxl ( 120021) for a complete re- 
view of unsolved problems in stellar pulsation physics. 

With two observed frequencies (Rodrig uez et alJ Tl992l) . 
the 6 Scuti star RV Arietis (GSC 01217-01057) belongs to 
the subclass known as high-amplitude 5 Scuti stars (HADS). 
Suc h characteristics w ere first attributed to RV Ari by 
iRodriguez et al.l dl992l) . after they analysed uvby/3 Stromgren 
photometric observations, and then identified the observed fre- 
quencies (Table [Q as the fundament al radia l mode and its first 
overtone. More recently, IPocs et all (2002) examined the am- 
plitude variation of both modes and found an additional fre- 
quency at fj, = 13.61 16 cd . However, the same authors claim 
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that more observations are needed in order to confirm this third 
frequency and we have not used it in the present work. 

Such stars have long been considered similar to double- 
mode Cepheids, i.e., evolved stars pulsating only in the first 
two radial modes. However, in the last years, low-amplitude, 
non-radial modes have been found together with the radial 
ones, making these stars more similar to the unevolved 6 Scuti 
stars, which are quite commonly distributed bet ween the zero - 
age and the termina l-age main se quence stars llPorettil E003). 
Nevertheless, Poretti et al ] ll2005t) have provided very recent 
evidence for the existence of pure radial double-mode pulsators 
within the HADS class and have stressed that they consitute a 
very homogeneous class of variable stars. 

IPoretti et all fcoOS) and lSuarez etalJ (f2006) studied RV Ari 
from a theoretical point of view and in the framework of the 
Petersen Diagrams. While the former considered it within the 
subset of classic Pop. I pure radial double-mode pulsators, the 
latter used it (along with other similar objects of the class) to 
study the impact of rotation on the determination of masses 
and metallicities as derived from the Petersen Diagrams. The 
independent information that such diagrams can provide turns 
out to be very helpful for modelling this kind of star. However, 
neither of these two works made a specific study of RV Ari. 

The present work aims at a comprehensive modelling of 
the HADS star RV Ari. To do so, some of most updated 
tools adapted for t his purpose are used: 1) the evolution- 
ary code CESAM dMorellll997l) and 2) t he pulsation codes 
GRACO JMova et alJ |2004 and FILOU jTran Minh & Lonl 
1 19951 ISuaredbOOaK The GRACO provides the non-adiabatic 
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Table 1. Observed mean periods of the three detected oscil- 
lation frequencies with their correspondin g mean amplitudes . 
The first two periods are reported by Rodriguez et al. 
and the third one is proposed by Po cs et alJ (J2002). The last 
column gives the frequency ratio f\ lfi = 2,i- 




Frequency P(d) v(c/d) v(pHz) A v (mag) fi/fi=2,i 



quantities related to pulsation and includ es the atmosphere- 
pulsation interaction described in Dupr et et all i200A The 
FILOU code includes the effects of rotation on adiabatic oscil- 
lations up to the second order in a perturbative theory. Using 
these tools we performed a massive numerical study of RV 
Ari to attempt to constrain physical and theoretical parame- 
ters with a precision never obtained for any other 6 Scuti star. 
Particularly, the following aspects are studied: 

• the effects of rotation on both the oscillation frequencies 
and the equilibrium models; 

• the use of the mixing-length theory for describing convec- 
tion in equilibrium models; 

• the calculation of non-adiabatic observables, needed for 
instability analysis, are also used for mode identification 
within the framework of multicolour photometry; 

• the use of updated Kurucz models for a detailed study of 
the atmosphere-pulsation interaction, 

• analysis of RV Ari in the framework of the Petersen dia- 
grams (PD), including the rotation effects. 

Such an exhaustive modelling intends not only to better con- 
strain the set of representative models of the star but also to in- 
troduce a comprehensive and self-consistent methodology for 
interpreting the oscillation spectra of 6 Scuti stars. 

The paper is structured as follows: equilibrium models and 
oscillation computations are described in Sect. |2] The funda- 
mental parameters of RV Ari and the location of the star in the 
HR diagram are given in Sect. [3] Then, a non-adiabatic analysis 
is performed in Sect. |4] including studies of mode-instability 
ranges and convective efficiency. Section[5]focuses on the prob- 
lem of mode identification, which is undertaken by consider- 
ing amplitude/phase diagrams in the framework of multicolour 
photometry, and examining the variation of colour indices as 
a function of the star position in the HR diagram. In Sect. [6] 
the effects of rotation are considered, and in Sect. [7J Petersen 
diagrams for this star are analysed. Finally, conclusions are re- 
ported in Sect. |8l 

2. The modelling 

In this section we describe the equilibrium models, the theo- 
retical oscillation computations and the linear stability analysis 
performed for the complete modelling of RV Ari. 
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Fig. 1. HR diagram containing the observational error box 
(shaded area) for RV Ari. Asterisks correspond to models rep- 
resentative of the star, labelled from (which corresponds to 
the observed mean values) to 4, covering the extreme regions 
of the error box (more details in the text). 



2.1. Equilibrium models 

The stell ar models wer e computed with the evolutionary code 
CESAM JMoretil997l) . A grid of around 2000 B-spline basis 
points was used in order to optimise the numerical precision of 
oscillation frequencies and eigenfunctions. 

The equation of state CEFF 

(Christensen-Dalsgaard & Daeppen] 1 1992b was used, for 
which the Co ulombian correction to the classical EFF 
( Egg leton et al.l fl973h was included. The pp chain, as well 
as the CNO cycle nuclear reactions, were considered and 
standard species from 'H to n O includ ed. Opacity table s 
were taken from the OPAL package Jlglesias & Rogers! 1996). 
complemen ted at low temperatures (T < 10 3 K) by the tables 
provided by Alexander & Ferguson ( 1994). In the atmosphere, 
the Eddington T(j) law (grey approximation) was considered 
for models used to compute rotating adiabat ic oscil l ations . 
For non-rotating, non-adiabatic calculations, a Kurucz (1998) 
model atmosphere was adopted. A metallicity of Z = 0.018 
was used, but allowed to vary when PD were analysed. 

Rotation effects on equilibrium models (pseudo-rotating 
models) were considered by modifying the equations 
( Kipp enhahn & Weigerl l 990 ) to include the spherically sym- 
metric contribution of the centrifugal acceleration by means of 
an effective gravity g e ff = g - 3l c (r), where g corresponds to 
the local gravity and J\ c {r) represents the radial component of 
the centrifugal acceleration. During evolution, models were as- 
sumed to rotate as a rigid body, and their total angular momen- 
tum was conserved. 

2.2. Oscillation computations & linear stability analysis 

Theoretical oscillation spectra were computed from the equi- 
librium models described in the previous section. For this pur- 
pose two oscillation codes were used: 
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Fig. 2. Predicted growth rates for modes in the range of the observed frequencies (left panels), and convective efficiency (defined 
in Sect. 14. 1 i as a function of the temperature (in logarithmic scale) along the star radius (right panels). Different line types 
correspond to different values of the mixing-length parameter a. The degree of modes are considered up to € = 3 and radial 
orders n are those required to cover the instability range. Vertical lines represent the observed frequencies f\ and /2. 
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Table 2. Observed amplitudes and phases (referred to y) and their corresponding uncertainties in the four filters of the Stromgren 
system. 



<Pu ~<Py <Pv~ <Py 


<pb ~ 




AJA y 


AJAy 


A„/A y 


fi 6.2 ± 2.0 1.8+1.7 


1.2 ± 


1.7 


1.086 ± 0.033 


1.394 ±0.037 


1.204 ±0.034 


f 2 3.8 ±6.7 -1.7 ±5.8 


-2.1 4 


5.7 


1.197 ±0.124 


1.492 ±0.139 


1.328 ±0.131 



1. The Granada Oscillation COde (GRACO, iMova et alJ 
120041) provided diagnostics on the instability and non- 
adiabatic observables required for fitting the multicolour 
photometric observations. In this code the n qn-adia batic 
pulsati on equations are solved mainly following lUnno etafl 
( 1989), including the non-adiabatic pulsation treatment of 
the atmosphere proposed by Dupr et et alJ (12002). Although 
including this atmosphere-pulsation interaction does not 
have any significant influence on the modal stability, it does 
allow theoretical predictions t o be compared w ith photo- 
metric colour observations (see Mova et al.l2004l for more 
details). The atmospheres were reconstructed at specific 
Rosseland optical depths (tr oss = 1) until the last photo- 
spheric edge of the star was reached. 

2. For rotating adiabatic oscillations computations, the oscil- 
lation code FILOU (Suarez 2002) provides adiabatic os- 
cillation modes corrected for the effects of rotation up to 
second order (centrifugal and Coriolis forces), including 
second-order effects of near degeneracy that are expected 
to be significant f or moderately high rotational velocities 
fSoiifi et alJ ll998). Indeed, two or more modes, close in 
frequency, are rendered degenerate by rotation under cer- 
tain conditions, corresponding to selection rules. In par- 
ticular, these rules select modes with the same azimuthal 
order m and degrees £ differing by or 2 JSoufi et all 
1998). If we consider two generic modes a\ = (n, I, m) and 
c*2 = {n',C ,m') under these conditions, near degeneracy 
occurs for \o- ai - cr ai \ < ctq, where o~ ai and a ai represent 
the eigenfrequency associated to modes a\ and ao, respec- 
tivel y, and crq represe nts the stellar rotational frequency 
(see Goupil et al.l2000l for more details). 



3. Locating RV Ari in the HR diagram 

The fundamental parameters of RV Ari (log T e g = 3 .859 +0.09, 
logg = 3.94+0.25, logL/L© = 1.11 ±0.25, [M/H]=0.01±0.20 
and vsim < 18 km s _1 ) were kindly provided by E. Rodriguez, 
who derived them from data reported in Rodriguez et al. (1992, 
2000; and references therein). This allows us to depict an error 
box in which our models can in principle be considered as rep- 
resentative of the star (see Fig.[0- On the other hand, rotational 
velocity, as well as the inclinati on angle, modifies the loca- 
tion of stars in the HR diagram. Mic hel et al J fl 999) propose a 
method of determining such effects on photome tric parameters. 
This method was further developed by Perez Herna ndez et al] 
i ll 9991) for 6 Scuti stars, showing that uncertainties around 
100 - 150 K in effective temperature and ~ O.lOdex in logg 
can be found for moderately rotating stars. That result was later 



confirmed by Sua rez et"al] (|2002). In the present case, consid- 
ering the absence of additional information on the inclination 
angle of the star, we considered an error box of ~ 150 K in T e s 
and ~ 0.25 dex in log g. Moreover, since RV Ari shows a rel- 
atively small projected rotational velocity, these uncertainties 
remain within the typical values for 6 Scuti stars. 

4. Non-adiabatic analysis 

Equilibrium non-rotating models were computed in the manner 
described in Sect. |2] In order to avoid redundant (and thereby 
unnecessary) model computations, the centre and the four cor- 
ners of the photometric error box, whose characteristics are 
summarised in Tabled were considered. Models are identified 
(ID) with labels ij with i = 0, ...,4 and j = 0, 3. The index 
i indicates the location in the HR diagram (0 j corresponds to 
central models represented in Fig.[0. For each 2, i.e., for each 
selected location, j indicates a different value of the convec- 
tive efficiency (MET) parameter a = cmlt used (see Sect.|3. 
This procedure allowed us to discard non-satisfactory models, 
thereby bounding the region in which models can be consid- 
ered as representative of RV Ari. 

4.1. Instability ranges & convective efficiency 

We basically attempted to significantly reduce the set of rep- 
resentative models by examining the region of the photomet- 
ric box where modes are predicted to be overstable. Predicted 
growth rates, 77 (positive for unstable modes), obtained for 
modes in the range of the observed frequencies, are depicted in 
Fig.|5](left panels) for the models considered. Since 77 is a dis- 
crete function of v, the instability ranges given in Table[3]were 
obtained by interpolating the value for which r\ changes its sign. 
According to this, almost all the selected models predict an un- 
stable observed frequency range. The main difference concerns 
the size of the region (in frequency) of unstable modes, while 
the modes at the upper boundary of the instability region may 
have very different frequencies. The total number of unstable 
modes decreases when the mass of the model increases. For 
the less massive ones, the predictions for a = 2 are quite dif- 
ferent to those obtained for a < 1. In particular, for the last 
value, large radial order n modes become stable. This ca n be 
expla ined by examining the convective efficiency defined JCoxl 
1980) as 



r = 



— (v rad - V) 



1/3 



(1) 



where V ra d and V represent the radiative and temperature gra- 
dients, respectively. Here, A basically represents the ratio be- 
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Table 3. Main characteristics of computed models representative of RV Ari. Labels ranging from to 4 indicate the location of 
the models in the HR diagram displayed in Fig.^ From left to right, M represents the stellar mass in solar masses M ; T e ff the 
effective temperature in K (on a logarithmic scale); g the surface gravity in cgs (on a logarithmic scale); X c the central hydrogen 
mass fraction; the age in Myr; a the mixing-length parameter; v m ; n -v max represents the frequency range (//Hz) of predicted 
unstable modes; Az/H p the size of the external convective zone normalised by the local pressure scaleheight measured at the 
bottom of the convection zone; Az/R t the fraction of stellar radius occupied by the external convective zone; Vf the fundamental 
radial mode (in //Hz); and finally, p represents the stellar mean density (in cgs). 



ID (ij) 


M 


To (]■ 


g 


X c 


Age 


or 


V 


v 


Az/H p 


Az/R, 


Vf 


p 


00 


1.75 


3.859 


3.94 


0.302 


1186 


2.0 


78 


380 


1.16 


0.33 


130.98 


0.192 


01 


1.75 


3.859 


3.94 


0.302 


1186 


1.5 


79 


336 


0.82 


0.18 


131.03 


0.192 


02 


1.75 


3.859 


3.94 


0.303 


1186 


1.0 


95 


355 


0.71 


0.14 


131.10 


0.192 


03 


1.75 


3.859 


3.94 


0.303 


1186 


0.5 


98 


361 


0.66 


0.13 


131.13 


0.192 


10 


2.10 


3.867 


3.70 


0.132 


860 


2.0 


71 


242 


0.81 


0.20 


81.82 


0.076 


11 


2.10 


3.868 


3.70 


0.132 


860 


1.5 


77 


250 


0.75 


0.18 


81.82 


0.076 


12 


2.10 


3.868 


3.70 


0.132 


860 


1.0 


78 


255 


0.74 


0.17 


81.82 


0.076 


13 


2.10 


3.868 


3.70 


0.132 


860 


0.5 


78 


267 


0.71 


0.16 


81.83 


0.076 


20 


2.01 


3.850 


3.69 


0.121 


985 


2.0 


63 


218 


1.03 


0.33 


81.32 


0.075 


21 


2.01 


3.850 


3.69 


0.120 


985 


1.5 


66 


216 


0.70 


0.18 


81.32 


0.075 


22 


2.01 


3.850 


3.69 


0.120 


985 


1.0 


66 


225 


0.70 


0.18 


81.32 


0.075 


23 


2.01 


3.850 


3.69 


0.120 


985 


0.5 


67 


224 


0.63 


0.15 


81.32 


0.075 


30 


1.57 


3.868 


4.19 


0.538 


850 


2.0 


117 


780 


1.95 


0.80 


208.59 


0.483 


31 


1.57 


3.868 


4.19 


0.538 


850 


1.5 


124 


491 


0.87 


0.16 


208.50 


0.483 


32 


1.57 


3.867 


4.19 


0.538 


850 


1.0 


127 


525 


0.74 


0.12 


208.59 


0.483 


33 


1.57 


3.867 


4.19 


0.538 


850 


0.5 


127 


551 


0.69 


0.11 


208.49 


0.483 


40 


1.49 


3.850 


4.19 


0.521 


1080 


2.0 


111 


770 


2.58 


2.71 


205.05 


0.485 


41 


1.49 


3.850 


4.19 


0.540 


985 


1.5 


119 


785 


1.83 


0.77 


211.31 


0.494 


42 


1.49 


3.849 


4.19 


0.542 


980 


1.0 


125 


383 


0.80 


0.15 


211.61 


0.494 


43 


1.49 


3.850 


4.19 


0.543 


975 


0.5 


125 


406 


0.67 


0.12 


212.00 


0.494 



As shown in Fig.|2](right panels), the convective efficiency 
obtained from models j and 2 j varies significantly with re- 
spect to the 4j models. On the other hand, models Aj present 
a very different convective efficiency for a — 2 compared to 
the other values. In general, the larger a, the larger the size of 
the convective zone, the more efficient the convection, and the 
more important the influence of convection in the pulsations. 
Intermediate stars in the temperature range of RV Ari, i.e., ap- 
proximately from 7380 K to 7080 K, show two thin external 
convective zones. High a values may give rise to extend the 
size of the convective zones artificially, which are then overes- 
timated. Such a scenario should be carefully considered since 
the limit of validity for MLT may be exceeded. 

In particular, insconsistencies appear when physical mag- 
nitudes are averaged over a distance aH p larger than the con- 
vective zone. In other words, the mixing-length theory remains 
valid when a < Az/H p (where Az is the extention of the con- 
vective zone, calculated using the Scwarzschild criterion, and 
H p the locally-defined pressure scaleheight, calculated on the 
basis of the convective zone). Otherwise, the convective eddies 
would run a longer distance than the convective zone itself. The 
values of Az/Hp, together with the fraction of the stellar radius 
occupied by the external convective zone, are listed in Tabled] 
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Fig. 3. Values of a used in models of Table [3] as a function of 
their corresponding Az/H p value. The straight line represents 
the limit of validity of MLT (see discussion in the text). 



tween the convective and radiative conductivities, and ao is a 
constant. 
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Fig. 4. Amplitude ratios and phase differences for models i3 (a - 0.5), i ranging from (top) to 4 (bottom), calculated for 
the observed frequencies f\ (left columns) and fi (right columns), respectively. Observed values are represented by asterisks, 
corresponding to the four Stromgren filters. 



for all the models. In order to be physically reliable, we can 
establish the equality a = Az/H p as the very limit of validity of 
MLT. This is represented by a straight line in Fig.|3j where a 
is displayed as a function of Az/H p . Therefore, it can be stated 
that models under the line are the only coherent ones with the 
MLT description, which corresponds to models computed with 
a = 0.5. Nevertheless, these results should be considered with 
care, since they may change depending on where H p is mea- 
sured in the convection zone. 

The ultra-precise asteroseismic observations of the forth- 
coming space mission COROT represent an important opportu- 
nity to study the sensitivity of these borders to a because of the 
larger number of frequencies expected to be detected at very 



low amplitude. This would provide a better understanding of 
the external convective zone when the MLT is used. 



5. Towards a mode identification 

Two of the three detected frequencies in RV Ari are identified 
as the fundamental radial mode and its first overtone. The only 
plausible argument to suggest such an identification comes 
from the ratio between frequencies, 0.772, which is typically 
found in Pop. I HADS. This possibility is examined in detailed 
in the framework of the Petersen diagrams in Sect. [7] In the 
following we use the available multicolour photometric obser- 
vations to provide additional information on the degree I of the 
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Fig. 5. Phase-amplitude diagrams showing the evolution of the fundamental radial mode (left panels) and first overtone (right 
panels) for the model 02 (a = 1). From top to bottom, the different panels depict the results for the Stromgren filters u, v, and b, 
referred to y. 



spherical harmonic associated to each observed pulsation fre- 
quency. 

5.1. Amplitude ratios vs. phase differences 

The linear approximation to non- radial flux variations of a 
pulsating star was fir st derived by Dziembowskil Jl977h and 
later reformula ted bvlBalon a & Sto bid dl979^) and bv IWatso"n1 
dl988h . Then, iGarrido et all dl990h show that "v" and "y" 
Stromgren filters can be used for distinguishing the degree I 
of the spherical harmonic associated to each observed pulsa- 
tion frequency. This technique, based on non-adiabatic calcu- 
lations, supplies amplitude ratios and phase differences for the 
different modes as a function of £ and is basically dependent on 
the modes and on the atmospheric physical conditions (more 



details in lMovall2003l iMova et alj f2004). In particular, pulsa- 
tion is highly non-adiabatic in stellar surface layers, in which 
thermal relaxation time is either of the same order or even 
lower than the pulsation period. Accurately determinating the 
eigenfunctions in these layers therefore requires using a non- 
adiabatic description that includes the entire atmosphere. This 
procedure then makes it possible to relate multicolour photo- 
metric observables with such eigenfunctions and therefore al- 
lows direct constraining of some unknown physical parameters 
through direct comparison with observations (see Table[5}- 

The computed models obtained in the previous section 
provide predictions not only about the region of the unstable 
modes but also on non-adiabatic quantities that are necessary 
for calculating phase differences and amplitude ratios for the 
different Stromgren colours (Sect. Theoretical predictions 
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Table 4. Main characteristics of computed models representative of RV Ari for the two rotational velocities considered. Labels i 
from to 4 indicate the location of the models in the HR diagram given in Fig.^ From left to right, Q. represents the rotational 
velocity in km s _1 , M the stellar mass in solar masses M , T e g the effective temperature in K (on a logarithmic scale), g the surface 
gravity in cgs (on a logarithmic scale), X c the central Hydrogen fraction, the age in Myr, Vf the frequency of the fundamental 
radial mode (in //Hz), and finally, vq represents the rotational frequency (in //Hz) of the model. 



ID M T eff g L X, Age v f v a 
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1.48 


0.135 
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22 


2.00 
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1.39 


0.126 
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1.57 
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42 


1.49 


3.850 
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202.46 
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Q, = 52.6 
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2.10 
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of both quantities are then compared with the observed ones 
obtained for f\ and fi- For the sake of clarity and brevity, di- 
agnostic diagrams are depicted in Fig. |4] only for the sets of 
models ii. The reference filter for amplitude ratios and phase 
differences is taken to be y. 

For each ij model, the theoretical frequencies used to build 
these diagrams are the closest to the observed ones, within 
an uncertainty of 5% (in frequency) approximately, which in- 
cludes effects coming from rotation, a and metallicity varia- 
tions around the observed value [M/H] = 0.01. Nevertheless, 
the uncertainties coming from rotation largely dominates. 

Let us first examine the results for the j models. Analysis 
of amplitude ratios and phase differences obtained for ( — 
modes reveals that the predictions for these modes are not sen- 
sitive to a. Accordingly, the best identification corresponds to 
the fundamental mode and first overtone. However, non-radial 
modes cannot be discarded. For instance, the observed fre- 
quency fx is found to be compatible with a I — 2 mode. Indeed, 
the smaller a, the closer the theoretical predictions. In contrast, 
the predictions for the radial modes are almost insensitive to 
a. Then, for this set of models, only non-radial modes could 
provide a possible discrimination of this parameter. 

For 1 j models, the predictions for f\ are quite different 
from those for fx. While none of the models reproduces the 
observed colour indices for the former, for the latter, four mod- 
els can be considered as representative. As for Oj, the smaller 
the value of a, the better the predictions. In particular, the best 
results are obtained for a = 0.5, with which fx is identified as 
I = 3. Multiple mode identification was obtained for fx, for 
which, nevertheless, a small set of solutions was obtained. 

Only one mode for the 2 j models is found to be compat- 
ible with observations for f\, which would identify it as the 
second overtone of the fundamental radial mode. In the case 
of fx, when considering low values for a, all modes are com- 



patible with observations within the error bars, so that a mul- 
tiple mode identification can be done. In the case of 3 j and 4 j 
models, we found again that the best predictions are found for 
models evolved with low a. Analysis of diagrams reveals that, 
for a < 1, /i and fx are identified only as t — 1 and £ = 3, 
respectively. 

The results obtained for the coolest models should be in- 
terpreted with care since the effective temperature of RV Ari 
is within the limit of validity of the frozen convection ap- 
proximation adopted in this work. At these temperatures, the 
interactio n between co nvec tion and puls ation described by 
Grigahce ne et alJ J2005I) and lDupret et alJ ll2005h may change 
the non-adiabatic eigenfunctions significantly. 

5.2. Colour index variations 

Studying the evolution of colour indices using the ampli- 
tude/phase diagrams can also be useful for constraining the 
evolutionary stage and eventually the metallicity of the mod- 
els. To do so, we examined the evolutionary tracks of Oj mod- 
els for which a mode identification can be reliably performed 
(see Sect. 12. 21 . In particular, the observed f\ and fx for 02 were 
identified as the fundamental radial mode and its first overtone. 
The evolution of the theoretical Stromgren indices u, v, and b 
for these two modes is depicted in Fig. [5] 

If we keep only the models lying inside the photometric er- 
ror box for each filter and identify f\ and fx as the fundamental 
radial mode and its first overtone within a +5% in frequency, 
it is possible to constrain the range of fundamental parameters 
given by the observations. Left and right panels indicate that 
both frequencies provide similar constraints on the age and ef- 
fective temperature of the models. In particular, the effective 
temperature and evolutionary stage are constrained to [3.849, 
3.860] dex and [1 190, 1270] Myr, respectively. Note that these 
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ranges are unusually narrow compared to the typical asteroseis- 
mic modelling of 6 Scuti stars. 

6. Effect of rotation 

The effect of the low-rotational velocity on the determination 
of fundamental parameters of RV Ari remains within the ob- 
servational typical uncertainties for 6 Scuti stars. Since no con- 
straints of its inclination angle have ever been provided, we 
will only consider two extreme values in the following for the 
rotational velocity of this star: 18kms~', which represents the 
minimum possible value (i.e., i = 90°) and 52.6 km s~', corre- 
sponding to an angle of i = 20°. These values reasonably cover 

the range of rotational velocities observed in HADS. 

I ML I I 

As shown by Dziembowski & Goode ( 1992), by Soufi et al. 

Jl998l) . and more recently bv |Suarez et alJ fcOOa) . the effect 
of rotation on adiabatic oscillation modes should not be ne- 
glected as regards, for instance, the effe ct of near degener - 
acy (see Sect. 12. 2> or differential rotation llSuarez et al. 2006). 
Indeed, near degeneracy not only affects the value of the os- 
cillation frequency, but also the corresponding angular compo- 
nents are somehow mixed up. As explained in Sect. 12.21 near- 
degenerate modes cannot be represented by a single spherical 
harmonic, but a linear combination of them is involved. One 
of the most important consequences of this physical effect con- 
cerns the identification of mode degrees using multicolour pho- 
tometry by means of the phase/amplitude diagrams discussed 
in the previous section. This was stud ied for the first time 
bv lDaszvhska-Daszkiewicz et alJ J2002h - who claimed that, for 
near-degenerate modes, such diagrams depend on the rotational 
velocity and on the inclination angle of the star. 

In order to take the aforementioned effects into account, we 
first built pseudo-rotating models representative of RV Ari in 
the manner described in Sect. |2] for the two rotational veloci- 
ties selected. The main characteristics of the selected pseudo- 
rotating models are listed in Table |4] Notice that the same la- 
bels are used to represent the different models for both the 18 
and the 52.6 km s _1 rotational velocities. This was conveniently 
chosen for the sake of clarity (the location of models evolved 
with a rotational velocity of 18 kms -1 and the equivalent with 
52.6 km s~' are quite close in the HR diagram). 

6.1. Near degeneracy and phase-amplitude diagrams 

Using the pseudo-rotating models described in Table|4]we com- 
pute the corresponding oscillation frequencies in the manner 
described in Sect. I2.2I In principle, the presence of non-radial 
modes (up to £ — 3) is not discarded and the investigation is 
extended to degrees up to £ < 3. Higher I values are thought 
to suffer from cancellation effects, thereby affecting its visi- 
bility. As expected, for models with a rotational velocity of 
18kms _1 , the effect of rotation on the oscillation frequencies 
(including the effect of near degeneracy) is found to be small 
as far as the position in the HR diagram, metallicity, or a value, 
is concerned. Thus for this rotational velocity, the results given 
in the previous sections remain valid. However, if the star is 
not observed from the equator, for example with i = 20° (i.e., 
Q = 52.6 km s~'), the effect of rotation cannot be neglected. 



Table 5. Coupled modes found for the z'2 models with Q = 
52.6kms~ I (see Table|4}. The mode identification is given by 
its radial order n, spherical degree {, and azimuthal order m. cr a 
and (Tb (in yuHz) represent the oscillation frequency. /3j is the 
contamination coefficient (see text for definition). 
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Since we are interested in studying the impact of the rotation ef- 
fects in the whole method and, more particularly, in the phase- 
amplitude diagrams, only models rotating at Q. = 52.6 km s~' 
are considered hereafter. 

A complete list of degenerate modes close to the observed 
frequencies is reported in Table |5] We calculate the ratio be- 
tween their horizontal components for each degenerate mode 
couple. Following the formalism given in Suar ez et alJ (12006. 
hereafter SGM) for two coupled modes j = (a, b), this relation 
is written in the form: 



(1) , (2) , SU)Q OOJ Q 
Ka + m aa + ~TT + 

2 swo 

where m^V and rrl?) (i, j = a, b) represent the elements of the 
first- and second-order matrix Al,/ given, respectively, by Eqs. 
(22) and (23) of SGM, obtained from the oscillation equation 
in the presence of near degeneracy (SGM, Eq. 20). The non- 
perturbed eigenfrequency of modes are included in the terms 
wo and 5uq as defined in SGM (Eqs. 16 and 17, respectively). 
Hereafter, fij quantities are called contamination coefficients. 
They roughly provide the ratio of the original horizontal com- 
ponent in the coupled mode. Analysis of the contamination co- 
efficients of modes listed in Table|5]reveals that, systematically, 
the lowest value of /3j is assigned to the lowest t value and 
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Fig. 6. Phase-amplitude diagrams for the Stromgren filter u (relative to y), for theoretical near-degenerate modes whose fre- 
quencies are close to the observed ones. The squares represent the location of modes obtained when rotation is not taken into 
account. 



vice-versa. Inspired by genetics, such modes can be called re- 
cessives, and those with the highest /3j value, dominants. 

As claimed in iDaszvriska-Daszkiewicz et"all ll20Q2h . the 
coupled-mode positions in the amplitude/phase diagrams be- 
come both aspect- and m-dependent. Following that work, we 
calculated amplitude/phase diagrams for modes of Table [5] 
varying the angle of inclination i from -180° to 180°, in steps 
of cos i = 0.02°. The results are depicted in Fig. [6] in which 
only the filter u (as compared to y) is used. For the sake of clar- 
ity and brevity, only recessive modes are considered, because 
dominant modes show similar results and do not provide addi- 
tional information. For comparison, the values of the amplitude 
ratios and phase differences obtained for the non-rotating case 
are also depicted. Note that the density of points varies along 
the curves. This means that the amplitude and phase predictions 
are found to be very close each other for a significant number of 
inclination angles. Points of curves corresponding to I — 2 and 
3 modes are systematically located close to i = ±180° regions. 
The concerned dominant modes present contamination coeffi- 
cients close to 1 ; therefore, their contribution to the horizontal 
component of the coupled mode is almost the maximum possi- 
ble. However, if the star is observed at an intermediate i angle, 
the situation becomes less dramatic, and it is possible to iden- 
tify the recessive modes correctly. In the present case, it means 



that the ( identification performed using amplitude/phase dia- 
grams without taking the rotation effects into account remains 
valid for ( — and 1 modes. 

7. Petersen diagrams for RV Ari 

The Petersen diagrams (PD) show the ratio between the fun- 
damental radial mode and the first overtone periods, Ili/o, as a 
function of the fundamental mode period, logflo, for a given 
stellar model. These well-known diagrams are particularly use- 
ful for constraining the mass and metallicity of double-mode 
radial pulsators like HADS. In the case of RV Ari, the ratio be- 
tween the observed frequencies fa and fa is 0.772, which lies 
in the range of values typically observed in Pop. I double-mode 
pulsators (see Sect.^. Our previous modelling points towards 
the possibility that fa and fa correspond to the radial fundamen- 
tal mode and its first overtone (for 0/ models). It is thus worth 
investigating the observed periods in the framework of the PD. 

In Fig. a typical PD for z'2 non-rotating evolutionary 
tracks is depicted. The shaded region indicates the typical 
IIi/o values found for Pop. I stars in the range of IIi/o = 
[0.772,0.776]. The five tracks were computed with the ob- 
served metallicity [M/H]=0.01 dex, for which the range of 
masses [1.49, 2.10] M Q was found. The mass for RV Ari is then 
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Fig. 7. Typical PD (IIq in d) containing different tracks that cor- 
respond to the evolutionary sequences of z'2 models, computed 
with a metallicity of [M/H]=0.01. The shaded area corresponds 
to typical values found for Pop. I stars. The filled symbol rep- 
resents the observed Ili/o for the double-mode high-amplitude 
5 Scuti star RV Ari. 



Fig. 8. PD (Flo in d) showing evolutionary sequences of the cen- 
tral non-rotating model 02, computed for 5 different metallic- 
ities around the observed [M/H]=0.01. As in Fig.0 the filled 
symbol represents the observed IT /o for RV Ari, and the shaded 
area corresponds to typical values found for Pop. I stars. 



rpd, n 



predicted between 1.75 and 2.01 M n , in agreement with the re- 
sults obtained bv IPoretti et alJ (J2005) when considering solar 
metallicity. 

Furthermore, PD diagrams can also be used to constrain 
the me t allicity of representative models. As shown bvlPetersenl 
Jl973h : IPetersen & Christensen-Dalsgaard \ 19961 1 19991) . the 
period ratios increase when decreasing the stellar initial metal 
content. This property is illustrated in Fig. |8j where tracks of 
the non-rotating model 02 are depicted for the five metallici- 
ties studied in the previous sections. Accordingly, each track 
corresponds to models with a different mass, from 1.62 M 
for [M/H]=-0.24dex to 1.86 M for [M/H]=0.26dex. From 
this figure, the best agreement between models and observa- 
tions is obtained for the track of [M/H]=-0.24dex, for which 
the closest model has a mass of 1.62 M , log Tes = 3.858, 
logL/Lo = 1.09, and an age of 1478 Myr. In principle, this 
would correspond to 3 j models (see Table |3 rather than j 
models. To fit the j models, it would thus be necessary to in- 
crease the value of the metallicity of the models (and therefore 
the mass) and/or m odifying their rotational velocity. Recently, 
ISuarez et alJ J2006I) showed that rotation can also modify the 
PD, and must be taken into account when performing accu- 
rate predictions for these parameters, e ven for slow rotators. 
Indeed, as shown in Sua rez et al l ll2006l) . rotation increases the 
IIi/o period ratios for a given mass and metallicity. In contrast, 
for a given mass (and rotational velocity), such period ratios in- 
crease when decreasing the stellar metal content. This is illus- 
trated in Fig. [5] where, for the j models computed with the ob- 
served metallicity, we compare the non-rotating track shown in 
Fig.Qwith two other pseudo-rotating models with Q = 18 and 
52.6kms~ 1 . As expected, when considering Q = 18kms _1 , 
models are predicted to be located in the PD very close (larger 
ITi/o ratios) to those in the absence of rotation. On the other 
hand, if we consider Q. = 52.6 km s~', models are predicted 
to show higher period ratios than the observed ones. Thus, we 
can argue that the best location would be somewhere in be- 
tween, that is, in the range [1.75, 2.01] M Q with a rotational 




Fig. 9. PD (Flo in d) showing evolutionary sequences of the cen- 
tral pseudo-rotating model 02, computed for the two rotational 
velocities considered in this work: 18 and 52.6 km s~'. A third 
track corresponds to the non-rotating case shown in Fig.0 



velocity around 35 km s 1 (which corresponds to i = 30° ap- 
proximately) and with a metallicity approximately in the range 
[-0.09,0.11]. 



8. Conclusions 

We have performed a comprehensive theoretical study of the 
HADS star RV Ari using some of the most modern tech- 
niques available for modelling intermediate-mass pulsating 
stars. Updated Kurucz models were included in order to take 
the atmosphere-pulsation interaction into account. Different 
approaches were followed to calculate the oscillation frequen- 
cies. Non-adiabatic observables were calculated in order to per- 
form instability analysis and mode identification within the 
framework of multicolour photometry. Moreover, adiabatic os- 
cillation spectra corrected for the effect of rotation up to the 
second order (in a perturbative theory) were also calculated. 
This allowed us to examine the effect of rotation on previous 
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studies, in particular, for multicolour diagnostic diagrams and 
for Petersen diagrams. 

Non-adiabatic calculations point toward models located at 
the high-luminosity regions of the observational error box. 
Detailed analysis reveals strong constraints on the efficiency of 
the convection. In particular, we showed that the mixing-length 
parameter, a, is constrained to values close to a = 0.5. 

Multicolour diagnostic diagrams are found to be compati- 
ble with the assumption that the observed frequencies f\ and 
f2 are radial. However, the error bars of the observed ampli- 
tudes and phases for the different Stromgren filters do not allow 
modes to be identified in an unmistakable way. Nevertheless, 
all the amplitude ratios variations are predicted to be compati- 
ble with the calculated theoretical models that fulfill the above 
constraints. The study of colour index variations has allowed 
us to significantly constrain the effective temperature and evo- 
lutionary stage of the models to [3.849, 3.860] dex and [1190, 
1270] Myr, respectively. This constitutes one of the most accu- 
rate modellings ever performed for a 6 Scuti star, and for RV 
Ari in particular. 

As expected, rotation was found to affect the ampli- 
tude/phase diagrams depending on the inclination angle i of the 
star for near-degenerate modes. However, analysis of the vari- 
ations of amplitude ratios and phase differences as a function 
of i (using the contamination coefficients) has revealed that the 
amplitude/phase diagrams are still useful for the most probable 
i configurations of the star. 

Petersen diagrams were studied including the effect of ro- 
tation. Frequency ratios were found to be compatible with 
the theoretical predictions obtained with other modelling tech- 
niques discussed previously. 

Interestingly, the coherence and completeness of the results 
obtained (see Table |6j constitute strong evidence of the capa- 
bilities and robustness of the method presented here. Indeed, 
the full capabilities of the method will be tested once more 
frequencies with lower amplitudes ha ve been detected b y the 
forthcoming space mission COROT 1 (Bagli n et al .120021) . 
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Table 6. Compilation of the most restrictive ranges of physical parameters obtained in the present work for RV Ari. The first 
line corresponds to the results obtained using non-rotating, non-adiabatic models, and the second line shows results from using 
rotating, adiabatic models. 





M 


log r eff 


log L/L Q 


logg 


Q 


Age 


[M/H] 




M 


dex 


dex 


dex 


kms~' 


Myr 


dex 


NRNAM 


1.70-2.00 


3.860 - 3.849 


1.13- 1.14 


3.94 - 3.89 





1190- 1270 


0.01 ± 0.20 


RAM 


1.75-2.01 


3.903 - 3.821 


1.38-1.15 


3.93 - 3.77 


18-53 


795 - 1430 


-0.09-0.11 



